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GEOMETRY  AND  ANALYSIS  OF  SOME  PARETO-EFFICIENT 
EMPIRICAL  PRODUCTION  POSSIBILITY  SETS 

by 

A.  Charnes 
W.  W.  Cooper 

Z.  M.  Huang 
J.  J.  Rousseau 

Abstract 

Efficient  facets  for  the  additive  model  of  DEA  (Data 
Envelopment  Analysis)  are  obtained  by  a  new  series  of  linear 
programming  models  which  are  used  to  show  how  rates  of  change  of 
outputs  with  respect  to  inputs,  both  analytically , and  compu¬ 
tationally,  may  be  determined  along  a  given  efficient  facet.  The 
differences  between  single  and  multiple  outpu;  cases  are  determined 
and  exemplified  via  a  "cone  direction"  development. 

Keywords:  Data  Envelopment  Analysis,  Efficiency  Facet  Tradeoffs, 
Production  Efficiency  Tradeoffs,  and  Tradeoff  Geometry 
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EFFICIENT  FACETS  AND  RATES  OF  CHANGE: 
GEOMETRY  AND  ANALYSIS  OF  SOME  PARETO-EFFICIENT 
EMPIRICAL  PRODUCTION  POSSIBILITY  SETS 

by 

A.  Charnes 
W.  W.  Cooper 
Z.M.  Huang 
J.  J.  Rousseau 


1.  Introduction 

A  variety  of  models  have  now  been  developed  within  the  methodology  of 
Data  Envelopment  Analysis  (DEA)  for  evaluating  the  (relative)  efficiency  of 
Decision  Making  Units  (DMUs)  which  use  multiple  inputs  to  produce  multiple 
outputs. 

Building  on  the  engineering  ratio  idea  of  a  single-input,  single-output 
efficiency  measure,  Charnes,  Cooper  and  Rhodes  (1978,  1981)  generalized 
this  notion  to  multiple-input,  multiple-output  situations  and  pushed  forward  on 
both  their  managerial  aspects  and,  on  the  dual  side  of  the  mathematical 
programs  involved,  to  more  classical  notions  of  Pareto  efficiency  or  optimality. 
Subsequent  extensions  and  elaborations  in  DEA  include  the  "multiplicative" 
models  of  Charnes,  Cooper,  Seiford  and  Stutz  (1982,  1983),  the  "modified 
ratio"  model  of  Banker,  Charnes  and  Cooper  (1984),  the  "additive"  model  of 
Charnes,  Cooper,  Golany,  Seiford  and  Stutz  (1985),  and  the  "extended 
additive"  model  of  Charnes,  Cooper,  Rousseau  and  Semple  (1987).  These 
works  provide  alternative  (but  related)  notions  of  efficiency  with  respect  to 
different  empirically-defined  production  possibility  or  reference  sets.  Section  2 
and  Figure  1  summarize  the  construction  of  these  different  reference  sets. 
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The  identification  of  these  models,  indeed  the  whole  of  DEA,  as  Charnes- 
Cooper  (1961 ,  ch.  IX)  tests  for  Pareto-efficiency  of  the  generators  of  the 
empirically  defined  production  possibility  set  was  made  by  Charnes,  Cooper 
and  Seiford  (1981)  and  elaborated  in  Charnes  et  al  (1985).  They  also  began 
investigation  of  (among  other  informatics  and  function  properties)  the  efficient 
facets  of  empirical  production  possibility  sets  as  a  step  toward  determining  the 
rates  of  change  of  outputs  with  changes  in  inputs  or  other  substitutional  rates 
along  those  facets.  These  rates  of  change  have  important  economic  and 
managerial  implications  in  trade-off  analysis  and  resource  allocation. 

In  the  present  paper  we  extend  and  develop  in  depth  the  insights  and 
beginnings  in  Charnes  et  al  (1985).  Section  3  provides  a  theoretical  basis  for 
identification  of  the  efficient  facets  of  the  empirical  production  possibility  set  that 
underlies  the  additive  model  of  DEA.  The  efficient  facets  may  be  obtained  by 
solving  a  new  series  of  linear  programming  problems,  one  for  each  DEA 
efficient  observed  input-output  point.  Section  4  then  shows  analytically  and 
computationally  how  rates  of  change  of  outputs  with  respect  to  inputs  can  be 
determined  along  a  given  efficient  facet.  These  rates  of  change  are  computed 
from  a  linearly  independent  (in  the  inputs)  subset  of  the  facet  points.  They  will 
be  different  along  different  facets.  In  the  single  output  case  we  always  obtain 
nonnegative  rates  of  change,  but  with  multiple  outputs  this  is  not  guaranteed. 
The  "cone  direction"  development  in  Section  5,  in  either  the  ouput  space  or  the 
input  space,  show  us  what  combinations  of  associated  substitutions  are  needed 
in  order  to  obtain  nonnegative  rates  of  change.  A  simple  two-output,  two-input, 
ten-DMU  example  is  carried  throughout  the  paper  to  illustrate  what  is  involved. 
Concluding  remarks  that  are  given  in  Section  6. 
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2.  Empirical  Production  Possibility  Sets 

Consider  the  (empirical)  points  (xj ,  yj ),  j 1 , . . .  n,  where  the  xj  are  (mx1 ) 
input  vectors  and  the  yj  are  (sx1)  output  vectors.  In  most  applications  they  will 
be  positive  or  nonnegative  vectors.  We  define  the  'empirical  production  set', 

Pe,  to  be  the  convex  hull  of  these  empirical  points,  that  is, 

n  n  n 

Pe  =  {  (x.  y):  X  =  X  Xj  Pj,  y  =  X  Xj  4j.  V  pj  >  0,  X  Pj  =  1 }  (21) 

i-i  i-i  i-i 

as  shown  in  Figure  1 . 


Figure  1. 


The  'empirical  production  possibility  set’  Qe  of  Charnes  et  al  (1985)  is 
defined  by  adding  to  Pe  all  points  with  inputs  in  Pe  and  outputs  not  greater  than 
some  output  in  Pe  that  is, 

Qe  s  { (x,  y) :  X  =  X,  y  ^  y  for  some  ( x,  y)  e  Pe  }  (2-2) 

Thus  Qe  =  Pe  u  A  in  Figure  1 . 
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I  The  Banker,  Charnes  and  Cooper  (1984)  production  possibility  set  adds 

I 

to  Qe  the  set 

{ (x,  y):  X  ^  X,  y  =  y  for  some  (x,  y)  e  Qe  } 
is  given  by  Qe  u  B  in  Figure  1 . 


The  production  possibility  sets  studied  by  Farrell  (1957),  Shephard 
(1970),  and  FSre  and  Lovell  (1978)  are  truncated  cones,  given  by  Qe  uBuC  in 
Figure  1. 


For  efficient  production  we  wish  to  maximize  on  outputs  while  minimizing 
on  inputs.  Thus  we  set 


Qk 


1  ^k^s 

k  =  s  +  i,  1  =  1 . m 


for  (X,  y)  e  Qe 


where  the  yk  and  Xj  are  the  and  components  of  y  and  x.  A  Pareto-efficient 
(minimum)  point  for  gi  (x,  y), . . . ,  gs+m  (x ,  y)  is  a  point  (x*,  y*)  e  Qe  such  that 
there  is  no  other  point  (x,  y)  €  Qe  for  which 


Sk  (X,  y)  ^  Sk  (X*.  y*) .  ks  1 . s  +  m  (2.3) 


with  at  least  one  strict  inequality.  Evidently,  the  Pareto-efficient  points  of  Qe  are 
those  of  Pe,  hence  we  can  restrict  attention  to  Pe- 

Charnes  and  Cooper  (1961)  showed  (for  general,  multiple  goals 
functions  gk(x,y) )  that  (x*,y*)  is  Pareto-efficient  if  and  only  if  (x*,y*)  is  an  optimal 
solution  to  the  reduced,  single  mathematical  (goal)  program 

8-t-m 

min  X  gk(x,y) 

k.1 


gk(x,y)  ^  gk(x*,y*),  k  =  i , . . . ,  s  +  m 


(2.4) 


(x,y)€PE 

The  constraint  inequalities  in  (2.4)  for  a  test  point  (x*.  y*)  may  be  written  as 
y^y*.  x^x*. 
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which  are  the  envelopment  constraints  of  DEA  for  an  observed  input  vector  x* 
and  corresponding  output  vector  y*. 

3.  Pareto-Efficient  Facets  of  the  Empirical  Production 
Possibiiity  Set 

As  shown  by  the  locus  of  points  a  b  c  d  in  Figure  2,  the  Pareto-efficient 
empirical  production  frontier  is  segmented  into  facets  of  efficient  observed  input- 
output  points.  The  rates  of  change  of  outputs  with  respect  to  inputs  along  these 
efficient  facets  have  important  economic  and  managerial  implicaitons  for  trade¬ 
off  analysis  and  resource  allocation.  Our  development  in  Sections  4  and  5 
shows  how  these  rates  of  change  (which  will  be  different  along  different  facets) 
can  be  derived  from  the  observed  input-output  points  that  lie  on  the  efficient 
facets. 


In  this  section  we  provide  a  theoretical  basis  for  determining  the  efficient 
facets  and  a  practical  method  for  identifying  the  observed  points  on  a  facet.  A 
small  numerical  example  illustrates  what  is  involved. 


Figure  2. 
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Let  F  be  a  facet  of  Pe  contained  in  the  hyperplane  -a^  x  +  =  0  such  that 

-a^x  +  p\^0  for  all  (x,y)  €  Pe  (3.1) 

as  shown  in  Figure  2. 

Theorem  3.1 :  Let  the  relative  interior  of  F  in  the  hyperplane  -a''’  x  +  =  0 

be  non-empty,  and  let  at  least  one  other  observed  point  be  outside  F.  Then  F  is 
an  efficient  facet  if  and  only  if  a  >  0,  p  >  0. 

Proof:  By  our  assumption  on  F,  int  (Pe)  ^  0. 

(=> );  Let  (x,  y)  be  a  relative  interior  point  of  the  facet  F  contained  in 
-a"^  X  +  p\  =  0.  Assume  to  the  contrary  that  there  exists  a|  ^  0  for  some 
i  (i  =  1 , . . .  m)  or  Pr  ^  0  for  some  r  (r  =  1 . s). 

(i)  3  tti  ^  0  for  some  i. 

Consider  the  new  point  (x>_,  y^J  s  (x  -  y)  where  ej  is  a  vector  of  zeros  with  a 
1  in  the  i^^^  position. 

Then  for  any  X,  >  0  we  have 

-a^  Xx  +  P^  yx  =  X  +  Xoj  +  P^  y  ^  -a^  x  -i-  p^  y  =  0. 

Hence,  for  small  X  >  0  we  have 

(xx.  Vx)  €  Pe 
and 

xx^  X.  yx=  y. 

with  strict  inequality  holding  for  xj.  This  means  that  (x,  y)  is  not  efficient,  which 
implies  the  facet  F  contained  in  -a'*’  x  +  p\  *  0  is  not  efficient,  a  contradiction. 

(ii)  3Pr^0  for  some  r 


Consider  the  new  point  (x^,  yx)  =  (x ,  y  +  Xe^)  where  er  is  a  vector  of  zeros  with  a 
1  in  the  r^h  position.  Then  for  any  X  >  0  we  have 
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Hence,  for  small  ^  >  0  we  have 

(xx-  Vk)  €  Pe 


and 


xx=  X,  yx^  y 

with  strict  inequality  holding  for  yr .  This  means  that  (x,  y)  is  not  efficient,  which 
implies  the  facet  F  contained  in  -a"^  x  +  =  0  is  not  efficient,  a  contradiction 

(<=) :  Assume  to  the  contrary  that  the  facet  F  contained  in  -a^  x  +  =  0  with 

a  >  0  and  p  >  0  is  not  an  efficient  facet.  That  is,  for  some  given  relative  interior 
point  (x,  y)  of  F,  there  exists  another  point  (x,  y)  e  Pe  such  that 

x^  X  and  y  ^y 

with  at  least  one  strict  inequality  holding.  But  this  implies 
-a^x  + p^y  >-a''’x  + p^y  =0, 
which  contradicts  (3.1). 


Q.E.D. 
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Now  suppose  that  the  empirical  "input-output"  point  (xo.Vo  )  is  an  efficient 
point.  Consider  the  following  linear  programming  problem. 

min  0 

-a’’’  Xo  -I-  p^Vo  =  0 

-a'rxi  +  p\j^0.  j  =  1 . n  (3  2) 

e  +  p^e  =  1 
a>  0,  p>  0 

Figure  3 

(-a3 .  p3 ) 


Lemma  3.1 :  Let  (a,  p,  0)  be  a  feasible  solution  of  (3.2),  and  let 

-a*^  x  -t-  p'  y  »  0i  with  a'^  e  +  p'  e  ■  1  and  ^  0,  p*  ^  0,  i  ■  1 . k,  be  all  the 

facet  hyperplanes  passing  through  the  efficient  point  (x©  ,  yo)  over  Pe  u  A  u  B 

k 

(see  Figure  3).  Then  there  exists  {Xj}  with  Xj  ^  o  and  S  ®  1  such  that 

i.1 
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(a,  p,  0)  =  S  Xi  (a\  p',0i,) 

Proof:  Let  A  ={(x,y) :  -a'^x  +  p‘  y  <  0j,  i  =  1 . K}.  Since  (a,  p,  0)  is  a 

feasible  solution  of  (3.2),  the  hyperplane  -oc'’’  x  +  p^y  =  0  is  a  support  pinae  at 
(xo  .  yo  )  such  that 

-a^  Xo  +  p^yo  =  0  (3.3) 

and 

-a’’’  X  +  p  ‘  y  <  0  for  all  (x,  y)  e  A  (3.4) 

k  j 

First  we  show  that  there  exist  {Xj}  with  Xj^  0  such  that  (a,  P,  0)  =  S  P  ^i). 

i  >  1 

that  is, 

(-a,  p,  -0)  =  S  P'  “Qi). 

i- 1 

Assume  to  the  contrary  that  the  following  system  (3.5)  has  no  solution, 
(-a,  p,-0)  =  S  P'  ~Qi).  (3-5) 

i-1 

Xi>0.  i  =  1 . K. 


By  Parkas'  theorem,  there  exist  (d,  y,  f)  such  that 
-a' d  +  p'  y-  f  0j  ^  0,  i  =  1 , . . .  k 


and 

-a^d  +  p^y-f  0  >0 

By  (3.7),  (d,  y,  f)  must  be  a  non-zero  vector. 


(3.6) 


(3.7) 


(i)f  =  0 


Consider  (x,  y)  =  (xq,  yo)  +  (d,  y).  By  (3.6)  we  have 
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-a'^  X  +  p'  y  =  (-a*^  Xq  +  p'  Vo)  +(-a'^  d  +  p'  y) 

=  0i  +  (-a'^  d  +  p'  Y ) 

^  dj  for  all  i 

That  is,  (x,  y)  e  A.  But  by  (3.7)  we  have 

-a'*'  X  +  p^  y  =  {-<xJ  Xq  +  p^  Yo)  +(-a'''  d  +  p^  y ) 

=  6  +  (-a^  d  +  P^  Y ) 

>0 

This  contradicts  (3.4). 

(ii)  f>0 

Let  (x,  y)  =  (^,  y).  Then  by  (3.5)  we  have 

•T  _  i'''  .T  i''' 

-a'  x  +  p  y  =  (--a'  d  +  p  Y)/^^0i  for  alii. 

That  is,  (x,  y)  e  A.  But  by  (3.7)  we  have 

-a^x  +  p^y  =  (-a^d  +  ^  y)/1>Q 
This  contradicts  (3.4). 


(iii)f<0 


Let  (x,  y)  =  (xo,  yo)  +  (Xq  Vo  -  j)-  Then  by  (3.6)  we  have 

-a'  X  +  p'  y  =  2  (-a'  Xq  +  p'  yo)  +{-<x'  d  +  p'  y )/  (-f) 
^  2  0j  -  Gj  =  ©j  for  all  i 
That  is,  (X,  y)  e  A.  But  by  (3.7)  we  have 

-a’’’  X  +  p^  y  =  2  (-a''’  Xq  +  p^  yo)  +  (-a’’’  d  +  p^  y )  /  (-f ) 

>2  0-0 
=  0 

This  contradicts  (3.4). 

Hence  there  exist  {A,i }  with  Xj  >  0  such  that 

(a.  P.  0)  =  t  Xi  (a*  p‘  0i). 

i-1 

T  T  •1’ 

Since  a'  e  +  p  e  =  1  and  a'  e  +  p  e  =  1 ,  1=1 . k, 

we  have 

j  k  ^  k  .T 

1=a^e  +  p  e=^Xia*e+]^XjPe 
i-1  i-1 


Q.E.D. 

Theorem  3.2:  Let  (a*,  p*,  0*)  with  a*  >  0  and  p*  >  0  be  an  optimal  basic 
solution  of  (3.2)  (e.g.,  using  an  extreme  point  method,  such  as  the  simplex 
method).  Then  -a*^x  +  p*^  =  0*  is  a  hyperplane  containing  one  efficient  facet 
passing  through  the  efficient  point  (x©,  y©)  of  Pe. 
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Proof:  By  Lemma  3.1  the  feasible  region  of  (3.1)  is  the  convex  hull  of  the 

i  ‘T 

{ (a',  p ,  ©i) },  where  -a'  x  +  p  y  =  6i,  i  =  1, ....  k  are  all  the  hyperplanes 
containing  all  facets  passing  through  (xq,  yo) . 

When  using  an  exteme  point  method,  the  optimal  basic  solution  to  (3.2)  must  be 
one  of  the  { (a',  p',  ©i) },  say  a*=  >  0,  p*  =  p^  >  0,  and  ©*  =  ©^ .  By 

Theorem  3.1 :  -a*^x  +  p*^  =  ©*  is  a  hyperplane  containing  one  efficient  facet 
passing  through  (xq,  yo)  • 


Q.E.D. 

Theorem  3.2  provides  only  a  sufficient  condition  for  an  efficient  facet. 
Since  we  cannot  guarantee  a*  >  0  and  p*  >  0,  in  practice  we  may  use  small 
numbers  e  >  0  and  employ  the  following  linear  programming  problem  to  derive 
the  hyperplane  containing  one  efficient  facet  passing  through  the  efficient  point 
(xo.  yo)  of  Pe. 


Min  © 

-a^Xo  +  P^yo  =  e  (3,8) 

-a^Xj  +  p^yj  ^  ©,j  =  1, . . . ,  n 
e  +  p^  e  =  1 
ee,  (fe  ee. 

At  an  optimal  solution  (a*,  p*.  ©*)  to  (3.8),  all  those  observed  (efficient)  points  j 
which  satisfy  their  respective  constraints  as  equalities  also  lie  on  the  efficient 
facet  contained  in  the  hyperplane  passing  throught  (xo,  yo)-  Such  points, 
together  with  (xq,  yo).  constitute  a  subset  (but  not  necessarily  all)  of  the  facet 
members.  Applying  (3.8)  to  other  members  of  the  facet  will  generally  reveal 
additional  points,  and  there  will  be  duplication,  overlapping  and  "resting”  of 
these  various  subsets  from  which  the  facet  may  be  identified  by  reduction. 

Thus,  by  applying  (3.8)  to  each  DEA-efficient  point  in  turn,  all  efficient  facets  and 
their  member  points  can  be  identified.  This  procedure  requires  little  additional 
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computational  effort,  since  moving  from  one  efficient  opint  to  the  next  involves 
changing  only  the  first  constraint  of  (3.8)  with  everything  else  unchanged. 

To  illustrate  the  above  procedures  consider  the  following  two-output,  two- 
input,  10-DMU  example  with  data  as  given  in  Table  1. 


Table  1. 


DMU 

yi 

Outputs 

72 

X1 

Inputs 

X2 

1 

2 

1 

9 

9 

2 

3 

1 

12 

8 

3 

2 

2 

7 

12 

4 

5 

3 

6 

10 

5 

4 

4 

10 

5 

6 

3 

3 

8 

10 

7 

6 

6 

12 

10 

8 

8 

2 

14 

6 

9 

1 

6 

12 

12 

10 

3 

5 

8 

8 

The  Pareto-efficient  points  are  DMUs  4,  5,  7,  8  and  10.  Applying  program 
(3.8)  to  each  efficient  point  in  turn  produces  the  results  given  in  Table  2.  By 
reduction  we  see  that  there  are  two  efficient  facets:  (DMU  5,  DMU  8)  and 
(DMU  4,  DMU  5,  DMU  7,  DMU  10).  The  facets,  for  both  the  outputs  and  the 
inputs,  are  depicted  in  Figure  4. 

Table  2 


Program  (3.8) 
Applied  to  DMU 

4 

5 

7 

8 

10 


DMUs  Revealed  as 
Being  in  the  Facet 
4.5 
4.5 

4,5.7.10 

5,8 

4,5,10 
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4.  Determining  Rates  of  Change 


We  begin  this  section  with  some  definitions  and  lemmas  important  in  our 
subsequent  development. 

Definition  4.1  Let  S  be  an  arbitrary  set  in  The  affine  manifold  spanned 

q  q 

by  S  is  given  by  M  (S)  =  {x:  x=  S  XiX',  x'  e\Z  =  1 ,  q  >  1} 


Definition  4.2  Let  S  be  an  arbitrary  set  in  E^^^.  The  linear  subspace  spanned 

q 

by  S  is  given  by  L  (S)  =  {x:  x=  ,Z  XjX',  x'  e  S,  >.|  e  e\  q  >  1} 

i»1 

Lemma  4.1  Let  S  be  an  arbitrary  set  in  E^^^.  Then  for  any  xeS  we  have 


M  (S)  =  X  +  L  (S-x) 


Figure  5 
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Definition  4.3  Let  S  be  an  arbitrary  set  in  E^.  The  convex  hull  of  S  is  given 
by 

q  q 

K  (S)  =  {x:  x=  Z  X-jX',  X'  6  S,  Xi  ^  0,  Z  X.i  =  1 ,  q  ^  1} 

i*1  i«1 

Lemma  4.2  Let  S  =  {x"* . x^}  and  Rank  (x2-x'' . xk-x^)  =  p-1 ,  and  w.  I. 

0.  g.  let  x2-xi, ....  xP-xi  be  linearly  independent.  Then 

M  (xi . xP)  =  M  (x^  . . . ,  xi<) 

Proof :  Since  L  (x^-x"* . x^x^)  =  L  (x2-x'' . xP-x^) ,  then  by  Lemma 

4.1  we  have 

M  (x^  . . . ,  xk)  =  x"*  +  L  (x^xi, ....  x'^x'i) 

=  x^  +  L  {x2~xf, . . . ,  XP“X1) 

=  M  (xi, . . .  xP) 

Q.E.D. 

Lemma  4.3  Let  x  =  ^  X  x'.  Then  x  is  a  relative  interior  point  of  Kfx"*, . . .  x^)  in 
M  {x\  . . .  xk). 

Proof:  X  is  a  relative  interior  point  of  K(x‘', . . .  xP)  relative  to  M(x‘', . . .  xP).  Since 
M  (xi, . , .  xP)  =  M  (xi, . . .  xk)  and 

K  (x'*, . . .  xP)  c  K  (x"*, . . .  xk),  we  have  that  x  is  a  relative  interior  point  of 
K(xi, . . .  xk)  in  M(x‘', . . .  xk) . 


Q.E.D. 


17 


W.I.o.g.,  let  (x'l,  yi) . (xk,  yk)  be  all  the  generators  of  a  given  efficient 

facet.  Let  Rank  (x^-x'' . xk-x"*)  =  Rank  (x2-xi, . . . ,  xP-xi)  =  p-1 ,  and  let 

X  =  i  X  x'.  By  Lemma  4.3,  x  is  a  relative  interior  point  of  (x'', . . .  xk).  Hence,  if 
P  i-1 

we  can  determine  the  rates  of  change  of  the  outputs  with  changes  in  the  inputs  at 
X,  we  also  have  them  at  any  point  of  that  facet. 

Let  x-xU  i  X  (xi-xi)  =  X  P 
P  i-2 

and 

y-y'  =  (/-y'>  =  Y  p 

P  i-2 


where 

X  =  (x2-x^  ....  xP-x^) 
Y=(y2-yi,....yP-yi) 
and  ^  . p^^‘ 

Then 

x’^{x-xi)  =  x'^XP 
=>(X'^X)-iX'^(x-xi)=  P 
=>  y-yi  =  YCX'^  X)-1  X'^ix -xi) 
=  W(x-xi) 

where 


W  =  Y(XTX)-iXT. 


1  8 


For  any  r,  we  then  have  =  rW  (x-x^) 


(4.1) 


where  rW  denotes  the  r^h  row  of  W. 

Now  allow  a  small  change  Ai^  in  the  particular  input  i©  ,  with  all  other  input 
values  unchanged. 


Let 


Xi  +  Ai  -  x.l  =  y  a,  (xj  -xl ) 

'0  'o  'o  J  '  Iq  'o' 


and 


Xj-xl  =  X  Ctj(xj-xl),  i;ti 


j-2 


X  (AJ  =  (xi . xJ  +Ai, . Xm)T. 


Then  we  have 


and 


X  (Aj^)  -x’  =X  ocj  (xi-x^)  =  X  a 
i-2 


y  (Aj )  -yi  =X  aj{yi-y^)  =  Y  a 
j-2 


where 


a  =  (02 . ap)T 


Then 


a  =(X^X)-'  X^  (x(A,.)-x') 


and 


y  (A,J-y' =  Y(X V'  X^  {x(A,J-x') 


=  W(x(AJ  -xi). 


Now, 


Vr (Aj j-Vr  _  (yr(\)-y?)-  (y-Vr) 

\  ^io 


rW(x  (Ajj  -x^)-(yryr) 

_rW(x-xVWriAi^-(yr-yr^) 

Aio 

=  Wh„  (by  (4.1)). 


Thus,  we  have 
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If  Rank  (x2-x‘' . x*<  -  x"! )  <  m,  then  for  some  Iq  x  +  Aei^  will  no  longer  lie 

on  the  facet  for  sufficiently  small  A  >  0.  Therefore,  we  need  to  project  ej^  to  the 

subspace  L(x2-x‘' . xp-x^)  and  then  obtain  the  rate  of  change  along  this 

projected  direction.  (See  Figure  5) 

Let  P  =  X  {XJ  X)-1  X  T,  so  that  P  is  a  projection  operator  from  the  space 
to  the  subspace  L  (x2-x'' . xp-x'')  . 


Let 

Pe- 

=  ttt; where  ®io  has  a  1  in  the  io'"  position  and  zeros  elsewhere.  Then 

II  PSio  II 

we  have 


dVr  (x+pd'°) 
dp 


p  =  0 


=  rWd'“,  r=  1. . . .  ,s; L=.1, 


m 


(4.3) 


Note  that  if  Rank  (x^-x^, ....  x^-x"*)  =  m,  then  P  is  the  (m  x  m)  identity  matrix. 


Continuing  with  our  illustrative  example  of  the  previous  section,  recall  the 
two  efficient  facets  were  (5,  8)  and  (4,  5,  7, 10).  We  shall  now  determine  the  rates 
of  change  for  these  facets. 


Efficient  Facet  (5,  8): 


so  that 


X®  -  X®  =(  4  )  and  y®  -  y®  =  (  4^  ). 


Since  x8  -  x5  0  it  is  linearly  independent. 
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Thus 

X=(4)and  Y  =  {  4^  )  and  (X^  X)-i  =  ^. 

Therefore 


W  =  Y  (X‘^X)-1X‘^= 


8  2 
-4  -1 


P  =  X  (X^X)-1X'^=  ^ 

'  ^  17  4  1 


d^- 


i-  =  -2- 


II  Pe^  II  VTTT  U 


^2^_Pe2 - L 


PejII  VITll/' 


The  rates  of  change  are  given  by 


dyi  (x+pd^) 
dp 


p_o=iWd^  *i^»0.97 


dyi  (x+pd^) 
dp 


dyg  (x+pd^) 
dp 


dyg  (x+pd^) 
dp 


p 


p 


P 


-0 


-0 


-0 


=  iWd^  = -4= -0.97 

fTz 


=  2Wd^  =--4=--0.49 
V272 


.2Wd^  =  --4= -0.49 

fn 


Note  that  the  last  two  rates  of  change  are  negative;  we  will  address  this  situation 
in  Section  5. 


Efficient  Facet  (4,  5,  7,  10): 


It  is  easily  verified  that  xS-x^  and  x^-x^  are  linearly  independent. 


Hence 


and 


1 

3 


(X'^X)-i  = 


25’  75 
l2  ^ 
75’  900 


Therefore 


W  =  Y(X'^X)-iX'^ 


fl  l\ 

6’  3 

1  1 
\2’  5/ 


and 


P  =  X(X^X)-lX^  =  l2,2 
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Thus  the  rates  of  change  in  this  facet  are  given  by 
"  6  ’  ~  3 

3xi  ”  2  ’  3x2  ”  5‘ 


5.  Cone  Directions  for  Non-negative  Rates  of  Change 

The  rate  of  change  of  output  yr  with  respect  to  input  Xj,  as  given  by  (4.3),  is 
not  guaranteed  to  be  non-negative.  When  negative  rates  of  change  are  given  by 
(4.3)  we  may  employ  the  following  development. 

Note  first  that  the  "projections"  of  the  outputs  of  y^ . y*^  aipng  the 

direction  h  ^  0  are  given  by  the  single  outputs  y^  =  hV.  •  •  • .  y'^  =  Now 
consider  the  following  two  DEA  problems  (5.1 )  and  (5.2).  They  have  the  same 
input  data,  but  the  outputs  of  (5.2)  are  the  projections  (along  the  direction  h)  of  the 
outputs  of  (5.1),  and  (xo,  yo)  is  one  observed  point  in  {(x\  yi), . . . ,  (x^,  yn) }. 

min  -e''’s+  -e^s” 

1  yi  Xj  -  s+  »  y° 

-ZxiXj-s+  =-x° 

2  Xj  =1 

Xj,  s+ s“  ^0 
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min  -h^s+  -e^s“ 

Z  hV  Xj  -  h  V  =  h  V 

-ZxiXj  -s"  =-xo  ^5  2) 

LXj  =1 


Xj,  s*,s~  >0 

If  h  >  0  then  from  (5.1)  and  (5.2)  we  have  that  (x^,  yo)  is  efficient  if  and  only 
if  (x°,  hTyo)  is  efficient,  that  is,  (xo,yo)  is  efficient. 

Now  let 

Yh  =  (f-y' . yP-y') 

=  (h^y2-y') . h'^(yP-yi) 

=  h^(y2-y< . yP-yl) 

=  hT  Y 
and 


WH  =  Yh(XV‘X^ 

=  hPY(X^X)-lx''’ 


shTW 


We  then  have 
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dy  (x+pd') 
dp 


p  =  0 


=  Wud'  =  h^Wd',  i  =  1. 


m 


(5.3) 


Thus,  if  in  (4.3)  there  exist  some  r  and  i  such  that  rWd'  <  0,  we  need  to  find 
a  direction  h  >  0  such  that 


rWdi>0,  i  =  1,...,m  (5.4) 

Case  (i):  If  (5.4)  has  at  least  one  positive  solution,  then  it  is  in  the  cone  of 
directions  which  contains  nonnegative  rates  of  change  on  this  efficient  facet. 

Case  (ii):  (5.4)  has  no  positive  solution  and  Rank  (x^-x^ . x'<-x'')  <  m. 

The  hyperplane  -a*  x  +  p  y  *  0  which  contains  the  efficient  facet  has  been 

determined  by  the  linear  programming  problem  (3.8).  Now  set  some  suitable 

input  vector  x*  such  that  x2-xi . xP-xi,  x*-xi  are  linearly  independent,  and 

determine  the  corresponding  efficient  output  vector  y*  by  solving 

,T  • 

-a  X*  +  p  y  =  0  .  We  thus  have  a  new  efficient  point  (x*,  y*)  and  an  extended 
efficient  facet.  The  corresponding  new  X*-  Y*,  W*  and  P*  are  given  by 

X*  “  ( x2-xi, ....  xP-xf,  x*-xi) 

Y*-(y2-yi . yP-yi.y-yi) 

W*  =  Y*(X*Y)-ix‘'^ 

p*  =  X*  (X  Yh 


and 
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_ 


Pei 


where  ej  is  a  vector  of  zeros  with  a  1  in  the  position. 


Then  by  (4.3)  we  have 


dVr  (x+pd') 
dp 


_  =rW^,  r=1 . s:i=1 . m  (5.5) 


p  =  0 


If  the  rates  of  change  given  by  (5.5)  are  still  negative,  the  above  procedure 
is  repeated,  further  extending  the  efficient  facet  one  point  at  a  time.  The  process 
stops  after  at  most  m+  1-p  iterations. 

Case  (iii):  (5.4)  has  no  positive  solution  and  Rank  ( x2-x^  . . . ,  xl^-xi)  =  m. 

In  this  case  we  can  determine  directions  in  the  inputs  such  that  the  rates  of 
change  of  outputs  with  respect  to  inputs  are  nonnegative. 


We  need  to  find  a  direction  d  which  lies  on  the  subspace 


L  ( x2-x'' . xk-x'')  such  that 


dyr  (x+pd) 
dp 


p.o-WdSO, 


for  all  r  =  1 , . . . ,  s. 


Since  d  lies  on  L  (  x2-x'', . . . ,  x^^xi),  d  has  the  representation 
d  =  Pz,  for  some  z  e  E"’. 

This  means  we  must  have 


r  WPz  ^0,  for  all  r  =  1 . s 
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WP2  >  0. 


Then  for  any  direction  d  in  the  cone  A  given  by 


A  =  {d:  d  =  Pz,  WPz  >  0} 


we  have  the  rate  of  change 


dyr  (x+pd) 


.=  rWd^0,  forall  r  =  1, . . . ,  s. 


Wd^O 


Recall  that  for  the  efficient  facet  (5,  8)  of  our  illustrative  example,  the  rate  of 
change  of  the  second  output  with  respect  to  each  input  was  negative.  To  derive 

nonnegative  rates  of  change  along  this  facet  we  need  to  determine  a  direction 
h  =  0  such  that 


hTWd'sO,  i  =  1,2. 


Thus  we  have 


dy(xtpd^) 

dp 


h'^Wd^ 


17^272 


(68  hi-34  hj)  S  0 


dy  (x+pd^) 
dp 


(34  hi-17h2)  >  0 
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p  =  0 


2 

17/17 


which  imply  that  such  a  direction  must  be  in  the  cone 

2  hi  ^  h2. 

We  conclude  this  section  by  highlighting  the  distinction  between  the 
single-output  and  the  multiple-output  cases.  Charnes  et  al  (1985)  showed  that  if 
a  Pareto-efficient  empirical  production  function  has  only  a  single  output,  then  it  is 

an  isotone  function.  Hence,  if  Rank  ( x^-x^ . x'<-x‘')  =  m,  we  can  always 

obtain  nonnegative  rates  of  change.  If  Rank  ( x^-x"! . <  m,  we  can 

extend  the  facet  by  the  procedure  given  for  case  (ii)  above,  and  thus  will  be 
guaranteed  nonnegative  rates  of  change.  In  contrast,  as  has  been  shown  in  this 
paper,  the  multiple  output  case  is  considerably  more  complex. 

6.  Conclusion 

The  present  paper  has  extended  the  existing  theory  of  Data  Envelopment 
Analysis  to  develop  what  rates  of  change  of  outputs  with  changes  in  inputs  can 
be  determined  on  the  Pareto-efficient  facets  of  an  empirically  defined  production 
possibility  set.  These  rates  of  change,  which  will  be  different  on  different  facets, 
are  important  for  effective  management  of  the  resources  (inputs)  employed  to 
obtain  desired  feasible  outputs. 

The  efficient  facets  can  be  obtained  by  solving  a  series  of  linear 
programming  problems,  one  for  each  Pareto-efficient  observed  input-output  point. 
It  is  shown  that  the  rates  of  change  can  then  be  computed  from  any  linearly 
independent  (in  the  inputs)  subset  of  the  facet's  points.  For  the  single  output  case 
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we  always  obtain  nonnegative  rates  of  change.  The  multiple  output  case  is  more 
complex,  and  nonnegativity  is  not  guaranteed.  However,  the  "cone  direction" 
development  of  Section  5,  in  the  output  space  or  in  the  input  space,  shows  in 
what  directions  change  must  go  to  obtain  nonnegative  rates  of  change.  A  simple 
example  was  developed  and  carried  throughout  the  paper  to  provide  a  clearer 
understanding  of  the  geometry  of  the  empirical  Pareto-efficient  functions  as  well 
as  to  clarify  the  steps  in  our  procedures. 
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